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The  Naval  Postgraduate  School  open  circuit  oscillating 


flow  wind  tunnel  was  used  to  study  the  blowing  requirements 
to  maintain  an  attached  turbulent  boundary  layer  in  an  oscil- 
lating freestream  flow  with  an  adverse  pressure  gradient. 
Boundary  layer  separation  was  visualized  through  the  use  of 
tvfts.  Freestream  flow  oscillation  frequency  was  found  to 
have  an  effect  on  the  flowing  required  to  maintain  attached 
flow.  This  frequency  dependence  exhibits  characteristics 
which  suggest  resonant  behavior. 
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I . INTRODUCTION 


The  study  of  viscous  unsteady  flows  have  become  increas- 
ingly  important  in  recent  years.  Typical  examples  of  un- 
steady flovjs  are  found  in  helicopter  blade  aerodynamics  and 
in  gas  turbine  engine  internal  flows.  The  development  of 
advanced  design  concepts  in  helicopters,  ships,  and  propul- 
sion systems  requires  a basic  understanding  of  the  nature  of 
unsteady  flows.  The  complexity  of  the  problem  does  not  lend 
itself  readily  to  either  analytic  or  experimental  study, 
even  when  limited  to  flow  with  laminar  boundary  layers. 
Furthermore,  most  of  the  problems  of  interest  are  characterized 
by  turbulence  and  turbulent  boundary  layers,  which  greatly 
increases  the  difficulty  of  the  study.  Regardless  of  the 
difficulties  involved  the  field  is  of  such  importance  as  to 
warrant  the  expenditure  of  time  and  effort  to  widen  the  un- 
derstanding of  its  characteristics. 

The  present  study  was  to  experimentally  examine  the  ef- 
fect of  different  levels  of  blowing  on  turbulent  boundary 
layer  separation  in  an  oscillating  freestream.  The  effects 
of  changes  in  frequency  and  magnitude  of  the  freestream 
oscillation  on  the  blowing  requirements  to  maintain  an  at- 
tached boundary  layer  were  studied. 
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II.  BACKGROUND 


A.  NON-STEADY  FLOWS 

1.  Analytic  Studies 

Previous  analytic  studies  have  been  confined  mostly 
to  flows  which  have  a laminar  boundary  layer.  The  earliest 
was  Stokes  [Ref.  1]  study  of  the  doubly  infinite  flat  plate 
oscillating  in  its  own  plane  in  a fluid  at  rest.  Rayleigh 
[Ref.  2]  considered  the  second  order  effects  of  this  prob- 
lem. Schlichting  [Ref.  3]  expanded  Rayleigh's  work  to  in- 
clude the  boundary  layer  assumptions.  All  of  these  were 
simple  unsteady  flows  without  a mean  flow  or  pressure  gra- 
dient. Lighthill  [Ref.  4]  later  treated  the  case  of  small 
magnitude,  low  frequency,  sinusoidal  oscillations  super- 
imposed on  a steady  flow  about  a cylinder  and  a semi-infinite 
flat  plate.  Using  small  perturbati-  n theory  and  restricting 
the  solution  to  first  order  terms  he  determined  that,  at  low 
frequencies  the  flow  is  essentially  quasi-steady.  The  un- 
steady flow  had  the  same  characteristics  for  any  instantan- 
eous magnitt  l.e  of  freestream  velocity  as  a steady  flow  of 
the  same  velocity.  As  frequency  is  increased  a limiting 
frequency  is  reached  beyond  which  the  boundary  layer  reacts 
as  it  would  without  a mean  flow.  Lin  [Ref.  5]  analyzed  the 
high  frequency  case  for  large  amplitude  oscillations.  He 
found  that  the  governing  equations  became  essentially  lin- 
eal’- at  high  frequencies  allowing  separation  of  oscillatory 
and  steady  components.  Nickerson  [Ref.  f>]  expanded  this 

1’ 


work  to  include  higher  order  perturbation  terms.  Barriol 
and  Lucius  [Ref.  7]  used  numerical  methods  to  obtain 
asymptotic  solutions  to  the  boundary  layer  equations  for 
oscillating  flow  on  a semi-infinite  flat  plate  with  no 
pressure  gradient.  Their  solution  agrees  with  those  of 
Lighthill  and  Nickerson. 

The  previously  mentioned  studies  ^.ave  all  dealt  with 
laminar1  boundary  layer  flows.  The  analytic  approach  to  tur- 
bulent boundary  layer  flow  has  apparently  been  less  fruitful. 
Karlssen  [Ref.  8]  separated  the  fluctuations  into  periodic 
and  random  components,  and  by  averaging  over  a complete 
period  obtained  equations  similar  to  the  steady  turbulent 
boundary  layer  equations,  but  with  an  additional  fluctuating 
term.  Other  analytic  works  LRef.  9,  10]  have  used  linearized 
solutions  to  solve  the  pressure  distribution  on  an  airfoil 
at  low  angle  of  attack,  but  are  applicable  to  only  a small 
number  of  cases. 

2 . Experimental  Studies 

A great  deal  more  work  has  been  done  experimentally 
than  analytically,  but  there  are  still  large  gaps  in  the 
field.  Nickerson  [Ref.  6]  partially  verified  both  his 
analysis  and  that  of  Lighthill.  Miller  [Ref.  11]  studied  • 
the  transition  phenomena  on  a flat  plate  in  oscillating  flow. 
He  was  able  to  confirm  some  of  the  previous  analytic  pre- 
dictions. In  addition  he  determined  the  transition  Reynolds 
number,  turbulent  Reynolds  number,  and  turbulent  jntermit- 
tency  factor  for  oscillating  Blasius  type  flow.  Despard 
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[Ref.  12]  investigated  the  separation  of  a laminar  -boundary 

I : " ' 

layer  in  oscillating  flow.  He  proposed  that  the  definition 
of  separation  for  a laminar  boundary  layer  in  an  oscillating 
freestream  be  the  initial  occurrence  of  zero  velocity  or 
reverse  flow  at  some,  point  in  the  velocity  profile  through- 
out the  entire  cycle  of  oscillation.  He  was  also  able  to 
make  some  prediction  about  the  behavior  o^  the  separation 
point . 

Morrissey  [Ref.  13]  investigated  the  effects  of 
large  amplitude  flow  oscillations  on  the  heat  transfer  from 
a flat  plate  with  a turbulent  boundary  layer.  Jacobs  [Ref. 
14]  studied  the  effect  of  oscillating  mean  freestream  on 
the  turbulent  intensity  distribution  in  a turbulent  boundary 
layer.  These  two  studies  indicate  that  unsteady  flows  ex- 
hibit no  significant  alteration  in  the  character  of  the  eddy 
diffusivity  or  the  turbulent  intensity  distribution  when 
compared  to  a steady  turbulent  flow.  Banning  [Ref.  15]  in- 
vestigated the  pressure  distribution  on  an  airfoil  in  a 
turbulent  oscillating  freestream.  Others  have  studied  the 
effects  of  oscillating  flows  with  compressibility  effects 
included . 

B.  BLOWING  TO  AVOID  SEPARATION 
1 . Experimental  Studies 

Boundary  layer  control  through  blowing  has  been  wi  ui 
us  almost  as  long  as  the  boundary  layer  concept  itself.  The 
idea  arose  initially  from  the  use  of  slots  as  a lift  augmen- 
tation. device,  which  was  suggested  by  Lachmann  of  Germany  in 
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1918  and  later  tried  by  him  and  Handley-Page  in  England. 

Betz  [Ref.  16]  theorized  that  the  effect  of  the  slots  was 
to  accelerate  the  boundary  layer.  Baumann  [Ref.  17]  was 
led  by  this  interpretation  to  replace  the  air  passing  through 
the  slot  with  air  ejected  from  the  wing  interior.  This  pro- 
duced roughly  the  same  effect  as  the  slots  with  the  added 
advantage  of  being  able  to  control  the  eAect  by  control  of 
the  pressure  inside  the  wing. 

Until  the  1940’s  little  more  was  done  with  blowing 
except  some  experiments  which  used  it  as  an  alternative  for 
ailerons.  This  was  done  by  blowing  out  a slot  over  a short 
wing  section  to  induce  high  lift.  By  differential  blowing 
on  the  wing  a rolling  moment  was  produced  similar  to  that 
of  an  aileron.  In  1942-43,  triggered,  no  doubt,  by  the  war, 
experiments  began  to  be  conducted  in  France , Germany  and  the 
United  States,  followed  shortly  by  Britain.  Probably  the 
most  important  contribution  of  this  era  was  nade  in  1948  by 
Poisson-Quinton  [Ref.  18].  Prior  to  this  tim<_  the  common 
parameter  to  measure  blowing  was 


where  Q.  is  the  volumetric  flow  through  the  slot,  U the 
D 00 

mean  freestream  velocity  and  S the  characteristic  area  of 
interest.  Schwier  [Ref.  19]  had  shown  that  narrow  slots 
were  more  efficient  than  wide  ones.  Poisson-Quinton  verified 
this  and  suggested  that  a more  appropriate  blowing  parameter 
was 
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( 1 1 — 2 ) 


m.V . 
C = -.U 
M q„S 


where  m.  is  the  mass  flow  of  blown  air,  V.  the  jet  velocity 
3 J 

at  the  exit  assuming  isentropic  expansion  from  the  total 
pressure  and  temperature  inside  the  jet  supply  duct,  and 

• • , i ’ . r 

the  frees'tr'ean  dynamic  pressure.  Subsequent  experimental 
work  supported  the  use  of  this  parameter  %o  predict  the 
ef fectiveness  of  blowing. 

Unfortunately  all  the  early  efforts  in  boundary 
layer  control  through  blowing  were  hampered  by  the  great 
difficulties  involved  in  practical  applications,  Major 
obstacles  were  a supply  of  high  pressure  air  and  the  weight 
penalty  associated  with  the  ducting.  However,  with  the 
advent  of  gas  turbine  engines  for  aircraft  propulsion  a 
good  supply  of  air  became  available.  W-i  th  this  development 
and  advances  in  metal  alloys  and  other  lightweight  materials 
of  high  strength,  blowing  has  become  much  more  attractive. 
Since  1950  the  use  of  blowing  has  proceeded  at  an  accelerat- 
ing pace.  Unfortunately  little  work  has  been  accomplished 
in  the  area  with  which  this  study  is  concerned.  The  majority 
of  the  work  has  been  in  the  use  of  blowing  for  lift  augmenta- 
tion by  delaying  stall  or  by  providing  circulation  control. 
All  of  these  experiments  present  results  in  terms  of  CA»CD» 

and  C and  their  variation  with  C . Also  the  values  of  C,, 
m y U 

required  for  lift  augmentation  are  much  greater  than  that 
required  to  keep  the  boundary  layer  from  separating  under 
the  influence  of  an  adverse  pressure  gradient.  As  a result 
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there  has  been  little  interest  in  studying  the  amount  of 
blowing  required  to  maintain  streamline  flow.  The  only  in- 
formation to  be  gleaned  from  this  work  was  that  higher  pres- 
sure gradients  require  higher  Cu  to  overcome  separation. 

2 . Analytic  Studies 

Theoretically,  most  of  the  interest  lies  in  increas- 
ing lift,  rather  than  detailed  study  of  blowing  requirements 
to  maintain  attachment.  There  are,  however,  a few  notable 
exceptions.  Carriere  and  Eichelbrenner  [Ref.  2?]  developed 
a method  for  calculating  the  conditions  for  flow  re-attach- 
ment by  a jet  discharging  against  adverse  pressure  gradient. 
Their  analysis,  unfortunately,  was  heavily  dependent  on  the 
availability  of  empirical  profiles  from  experimentation  and 
several  arbitrary  assumptions.  Within  these  limits  .it  did, 
however,  provide  solution  to  cases  with  weak  pressure  gra- 
dients. Kozlous  and  Zyanyak  [Ref.  24]  have  recently  anal- 
yzed a laminar  boundary  layer  in  an  unsteady  incoming  motion 
around  a body  of  arbitrary  shape  with  either  suction  or 
blowing.  They  used  a six  degree  polynomial  to  describe  the 
boundary  layer  under  these  conditions.  Their  result  was  an 
integral  equation  which  could  be  integrated  directly,  or 
reduced  to  a quadrature,  depending  on  the  velocity  profile. 
They  were  able  to  verify  their  analysis  for  a symmetric 
wing  in  a start-up  situation.  This  can  be  extended  to  an 
oscillating  flow  by  varying  the  description  of  the  incoming 
motion.  Although  derived  for  lift  and  drag  predictions  in 
laminar  flows  with  blowing,  this  analysis  shows  great  promise 
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for  being  extended  to  turbulent  boundary  layer  blowing  re- 
quirement predictions.  The  greatest  difficulty  with  this 
extension,  or  any  prediction  of  blowing  in  turbulent 
regimes , is  that  much  of  the  energy  imparted  to  the  boun- 
dary layer  is  by  turbulent  mixing  which  defies  accurate 
modeling . 

C.  APPLICATIONS 

Recently  there  has  been  work  accomplished  toward  the 
application  of  boundary  layer  control  in  an  oscillating 
flow.  Englar  and  Williams  [Ref.  12,  23,  24]  have  applied 
boundary  layer  energization  through  blowing  to  augment  the 
lift  of  a submarine  stern  plane  and  a symmetric  airfoil  at 
angle  of  attack.  They  have  also  applied  tangential  blowing 
to  the  blades  of  a helicopter  rotor  system  (the  circulation 
control  rotor)  with  encouraging  results. 

There  has  not,  however,  been  any  detailed  studies  of 
the  actual  flow  over  the  surface  and  its  behavior  under  the 
conditions  of  high  angle  of  attack,  large  pressure  gradient, 
turbulence,  oscillating  freestream,  or  blowing  to  energize 
the  boundary  layer.  The  prediction  of  blowing  requirements 
under  these  conditions  will  for  the  foreseeable  future, 
rely  heavily  on  empirical  methods  based  upon  experimental 
results . The  present  study  is  an  attempt  to  provide  some 
of  this  experimental  data. 


III.  DESCRIPTION  OF  APPARATUS 


A.  OSCILLATING  FLOW  WIND  TUNNEL 
1 . General  Description 

The  experimental  work  was  conducted  in  the  lowspeed, 
oscillating  flow  wind  tunnel  located  in  the  Aeronautics 
Laboratories  of  the  Naval  Postgraduate  School.  This  wind 
tunnel  is  of  open  circuit  design*  with  a 24-inch  square  by 
223-inch  long  test  section.  A plan  view  of  the  tunnel  is 
presented  in  Figure  III.l.  The  tunnel  inlet  is  eight  feet 
square,  resulting  in  a 16:1  contraction  ratio.  Three  high 
solidity  screens  located  in  the  inlet  section  just  upstream 
of  the  nozule  produce  measured  freestream  turbulence  in- 
tensities of  0.261  to  0.413  percent  for  the  velocities 
encountered  in  the  present  work. 

The  wind  tunnel  drive  consists  of  two  Joy  Axivane 
Fans  in  series,  each  of  which  has  an  internal,  100  horse- 
power, direct  connected,  1750  rpm  motor.  The  fan  blades 
are  internally  adjustable  through  a pitch  range  of  25  to 
55  degrees,  providing  a wide  operating  range  of  test  section 
velocities.  Two  sets  of  variable  inlet  vanes,  located  im- 
mediately upstream  of  each  fan,  are  externally  operated  to 
provide  control  of  test  section  velocity.  These  vanes  are 
of  multileaf  design,  and  preswirl  the  air  in  the  direction 
of  fan  rotation  to  reduce  fan  capacity.  The  total  range  of 
tunnel  velocity  is  from  10  to  250  ft.  per  second. 
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2 . Rotating.  Shutter  Valve 


Two  fundamental  methods  of  creating  an  oscillating 
flow  environment  have  been  employed  in  the  past.  Nickerson 
[Ref.  6]  introduced  oscillations  by  oscillating  the  mode] 
in  a steady  flow  environment.  This  method  severely  re- 
stricts the  range  of  attainable  frequencies  because  of 
mechanical  complications,  and  also  introduces  measurement 
difficulties.  The  other  method  is  to  actually  oscillate  the 
flow  over  a stationary  model.  Hill  [Ref.  25]  used  a sliding 
abutter  to  impose  oscillations  on  the  freestream  but  was  re- 
stricted by  mechanical  limitations  to  low  frequencies. 

The  most  successful  method  of  obtaining  an  oscillat- 
ing flow  with  large  ranges  of  frequency  and  amplitude  was 
that  employed  by  Karlssen  [Ref.  8]  and  later  by  Miller  [Ref. 
11]  in  his  investigation  of  transition  phenomena.  A rotating 
shutter  valve,  immediately  downstream  of  the  test  section, 
is  employed  to  superimpose  a periodic  variation  of  velocity 
on  the  mean  flow.  The  method  used  in  the  present  investiga- 
tion is  identical  to  that  employed  by  Miller.  The  shutter 
valve  consists  of  four  horizontal  steel  shafts  equally  spaced 
across  the  test  section.  The  shafts  are  slotted  to  accommo- 
date flat  blades  of  various  widths , forming  a set  of  four  • 
butterfly  valves  spanning  the  test  section.  Figure  III.  2 is 
a photograph  of  the  shutter  valve.  Each  blade  is  driven  from 
its  immediate  neighbor  by  means  of  a timing  belt  and  pulley 
arrangement.  The  bottom  shaft  is  driven  by  a five  horse- 
power, variable  speed,  electric  motor,  through  a timing  b<_-lL 
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and  pulley.  An  intermediate  shaft  between  the  motor  and 
shutter  valve  permits  a wide  variety  of  pulley  ratios.  This 
drive  arrangement  provides  a frequency  range  from  0.4  cycles 
per  second  to  the  first  critical  frequency  of  933  cycles 
per  second.  The  electric  motor  presently  in  use,  however, 
restricts  the  oscillation  frequency  to  a maximum  of  240 
cycles  per  second.  The  amplitude  of  oscillation  is  controlled 
by  blade  width.  Test  section  closure  may  be  varied  from  25 
to  100  percent.  The  resulting  amplitude  of  oscillation  of 
test  section  velocity  is  a function  of  frequency,  mean 
velocity,  and  pressure  gradient.  In  this  investigation 
blades  producing  83.3  percent  closure  were  used  resulting  in 
a perturbation  range  from  5 to  25  percent  of  the  local  mean 
freestream  velocity.  Figure  III. 3 is  a picture  of  the  shut- 
ter valve  drive  arrangement. 

3 . Test  Section 

The  wind  tunnel  test  section  is  shown  in  Figure  III. 4. 
Continuous  pieces  of  two-inch  thick  aluminum,  24  inches  wide 
and  223  inches  long,  form  the  upper  and  lower  test  section 
walls.  Each  of  the  sidewalls  consists  of  three  two-inch 
thick  panels  of  stress  relieved  lucite.  For  this  investiga- 
tion the  central  sidewall  panel  on  the  opposite  side  of  the 
tunnel  from  the  control  ocnscle  was  replaced  wi.th  two-inch 
thick  plywood  to  facilitate  the  mounting  oi  instrumentation 
and  access  to  the  model  plenum.  The  Lucite  panels  on  the 
console  side  of  the  test  section  are  hinged  and  may  be 
raised  hydraulically,  providing  access  to  the  test  section. 
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Figure  III. 5 shows  the  test  section  with  the  door  open. 
The  heavy  construction  of  the  test  section  is  intended  to 
minimize  deflections  induced  by  rapid  changes  in  static 


pressure . 

Previous  test  on  the  tent  section  velocity  profile 
have  shown  the  velocity  variation  is  less  than  one  percent 
of  mean  to  within  three  inches  of  the  wall  [Ref.  12]. 

Figure  III. 6 is  a photographic  view  of  the  test  section  and 
control  console  with  the  model  in  place. 


B.  MODEL 

The  model  used  in  this  investigation  consisted  of  two 
plugs  placed  in  the  test  section  with  a flat  plate  halfway 
between  them.  Each  plug  had  a smoothly  curved  leading  edge, 
followed  by  42-inch  straight  section,  then  a diverging  sec- 
tion as  shown  in  Fig.  III. 7.  The  maximum  thickness  of  each 
plug  was  six  inches  and  the  plate  was  h inch  leaving  5-3/4 
inches  between  each  plug  and  the  plate.  The  diverging  sec- 
tion departs  from  the  horizontal  by  27°.  The  slant  length 
was  13  inches  and  the  characteristic  area  chosen  for  calcu- 
lation of  blowing  coefficient  was  the  plane  area  of  the 
diverging  section  which  was  305.5  square  inches. 

Each  plug  contained  a 1700  cubic  inch  plenum  chamber. 

The  entrance  to  the  plenum  was  a two-inch  diameter  hole  to 
which  the  blowing  air  supply  hose  was  connected.  The  blow- 
ing slots  were  23.5  inches  wide.  The  upper  slot  height  was 
.04],  inch  and  the  lower  was  .055  inch . The  cliff  e.  re  nee  in 
slot  heights  provided  information  on  the  blowing  coefficient 


requirements  as  a function  of  slot  height.  Because  of  its 
smaller  slot  height  the  upper  jet  wuald  need  less  mass  flow 
but  higher  jet  velocity  to  achieve  the  same  blowing  coeffi- 
cient as  the  lower  slot.  Figure  III. 8 shows  the  plenum  and 
slot  configuration.  Blowing  air  was  provided  by  a Carrier 
three-stage  centrifugal  compressor,  driven  by  a 300  horse- 
power General  Electric  induction  motor.  It  was  capable  of 
supplying  1900  cubic  feet  per  minute  at  29.5  psia.  Each 
plenum  was  supplied  and  metered  independently  by  a gate  valve 
in  the  three-inch  supply  pipe.  Each  pipe  had  an  orifice 
plate  with  a 1.8  inch  hole  diameter  for  use  in  mass  flow 
measurements . 

The  flat  plate,  constructed  of  a slab  of  phenolic  material 
24-inches  wide,  60-inches  long,  and  h inch  thick  had  a rounded 
leading  edge  in  order  to  avoid  leading  edge  separation.  The 
last  five  inches  was  tapered  and  hinged.  For  this  study  the 
hinged  portion  remained  at  zero  angle  of  incidence  with  re- 
spect to  the  flat  plate.  The  plate  was  mounted  at  the  center 
line  of  the  test  section  and  ran  from  a point  7.5-inches  from 
the  leading  edge  of  the  plugs  to  a point  five-inches  back  to 
the  end  of  the  plugs  as  shown  in  Fig.  III. 7. 

Figure  III. 9 shows  a schematic  of  the  test  setup. 


C.  INSTRUMENTATION  AND  CALIBRATION 
1.  Freestream  Sensors 
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freestream  dynamic  pressure.  A hot  wire  probe  29.5  inches 
aft  of  the  model  leading,  edge  and  midway  between  the  plate 
and  the  lower  model  section  was  used  to  determine  magnitude 
of  velocity  perturbations  during  oscillations.  The  probe 
was  connected  to  a locally  manufactured,  single,  channel, 
hot  wire  anemometer  which  was  connected  to  a Tektronix  555 
Dual-beam  oscilloscope  for  display  and  read  out.  Previous 
experimentation  [Ref.  12]  has  indicated  that  the  hot  wire 
circuitry  is  linear  with  freestream  velocity.  A short 
calibration  run  was  made  which  verified  this. 

The  frequency  of  oscillation  was  obtained  from  a 
magnetic  pickup  mounted  out-board  of  the  uppermost  shutter 
blade  shaft,  as  seen  in  Fig.  III.  2.  The  output  frequency 
was  read  on  a Dynascience  digital  counter. 

The  pressure  distribution  in  the  tunnel  was  measured 
by  a series  of  flush  pressure  ports  one  inch  apart  in  the 
flat  plate.  They  ran  from  a point  five  inches  in  front  of 
the  diverging  section  of  the  model  to  2.3  inches  from  the 
end  of  the  model.  The  ports  were  connected  to  a 40  tube, 

250  centimeter  upright  manometer  board.  The  median  reading 
of  the  oscillations  was  the  recorded  pressure  value. 

2 . Blowing  Coefficient 

The  mass  flow  of  blowing  air  was  measured  by  a pair 
of  orifice  plates  in  the  three  inch  supply  lines.  The  ori- 
fice plates  were  calibrated  in  situ.  This  was  accomplished 
by  connecting  the  blower  supply  lines  to  a 76.5  cubic  feet 
par  minute  rotameter  manufactured  by  Fischer  and  Porter 
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Company.  The  pressure  drop  across  the  plates  was  measured 
by, a 36  inch  U-tube  water  manometer.  The  results  of  the 
calibration  were  plotted  against  the  ASME  values  as  shown 
in  Fig.  III. 10.  Th' . differences  were  so  small  as  to  be 
negligible,  therefore  it  was  deemed  reasonable  to  extrapolate 
the  curve  along  the  same  form  as  the  ASME  curve  to  cover 
values  not  covered  by  the  rotameter.  Figure  I'll.  10  was  then 
used  to  find  the  volume  flow  of  air  for  any  value  of  pres- 
sure drop. 

The  pressure  downstream  of  the  orifice  plate  was 
measured  with  a 36  inch  U-tube  mercurcy  manometer  vented  to 
the  atmosphere.  Air  temperature  remained  essentially  con- 
stant from  the  compressor  to  the  plenums.  With  the  volume 
flow,  pressure  and  temperature,  the  mass  flow  could  easily 
be  calculated. 

The  blowing  jet  velocity  for  calculation  of  has, 
historically,  been  defined  as  the  velocity  at  the  exit  as- 
suming isentropic  from  some  pressure  and  temperature  inside 
the  supply  duct.  These  were  measured  with  a special  purpose 
Kiel  temperature  probe.  The  pressure  output  was  connected 
to  the  same  upright  manometer  board  as  the  pressure  distri- 
bution ports.  The  thermocouple  for  temperature  measurement 
was  connected  to  a Leeds  and  Northrup  portable  potentiometer. 

3 . Visualization  of  Separation 

Separation  lines  were  visualized  with  soft  twine 
tufts  placed  along  the  diverging  planes.  The  lines  of  tufts 
spanned  the  test  section  to  within  throe  inches  of  each 
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sidewall  and  were  spaced  at  one  inch  intervals  across  the 
test  section  and  two  inch  intervals  down  +>'-  plane  as  seen 
in  Fig.  III. 11. 


Photograph  of  Rotating  Shutter  Valve. 


Schematic  of  Test  Setup. 


,«M**~*»rr>*t*  .«  -■  if i«r «««:■  •*>  vt«  *>j|i^w*ir^  : 


Pressure  Drop  Across  Orifice  vs.  Volume 


Figure  III. 11.  Photograph  of  Diverging  Planes  from  Downstream. 
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IV.  TEST  PROGRAM 


A.  EXPERIMENTAL  PROCEDURE 

The  basic  test  program  was  to  find,  for  some  value  of 
freestream  dynamic  pressure,  the  value  of  blowing  coefficient 
that  was  just  sufficient  to  maintain  an  attached  boundary 
layer,  then  to  study  the  effects  of  varying  the  frequency 
of  oscillation  on  the  blowing  coefficient  required. 

The  first  problem  was  to  select  values  of  mean  free- 
stream dynamic  pressure  for  the  study.  The  dynamic  pres- 
sures chosen  were  those  which  indicated  5,  7.5,  and  10  cen- 
timeters on  the  micro-manometex  This  decision  was  based 
on  three  criteria: 

1.  Below  5 centimeters  the  atmospheric  wind  conditions 
at  the  inlet  began  to  affect  the  test  section  conditions. 

2.  Above  10  centimeters  the  fans  could  not  maintain  the 
pressure  at  higher  oscillation  frequencies.  This  could  be 
remedied  by  opening  the  fan  housing  and  altering  the  vane 
direction  but  this  would  have  vastly  increased  the  time  re- 
quired for  runs  and  would  have  had  an  effect  on  test  section 
conditions . 

3.  These  three  settings  were  easily  read  and  maintained 
during  tunnel  operations. 

These  three  dynamic  pressures  gave  values  of  mean  freestream 
velocity  of  about  96,  118,  and  138  feet  per  second,  respec- 
tively, although  this  var'ied  slightly  with  oscillation 
frequency. 
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The  initial  tests  were  dictated  by  the  wind  tunnel  set- 
up from  previous  work.  The  set  of  belts  and  pulleys  con- 
necting the  electric  drive  motor  to  the  shutter  valve,  at 
the  beginning,  corresponded  to  a frequency  range  of  2.48  to 
14.4  cycles  per  second.  Then  on  the  basis  of  the  previous 
test  results  the  range  was  expanded  outward  to  cover  fre- 
quencies of  1,6  to  16  cycles  per  second.  (See  Results  for 
detailed  explanation.) 

The  detailed  test  procedure  was  as  follows: 

1.  Start  tunnel  and  stabilize  at  desired  dynamic 
pressure . 

2.  Adjust  blowing  to  barely  maintain  attached  flow 
down  the  diverging  plane. 

3.  Take  data  for  pressure  distribution,  blower  plenum 
conditions  and  mass  flow  of  blown  air. 

4.  Start  shutter  valve  and  set  to  the  desired  frequency. 
Adjust  dynamic  pressure  as  necessary. 

5.  Repeat  step  2. 

6.  Take  data  for  magnitude  of  oscillation. 

7.  Repeat  step  3. 

8.  Adjust  shutter  valve  for  new  frequency  and  repeat 
steps  5,  6,  and  7 in  order. 

B.  DATA  REDUCTION 

The  data  reduction  was  accomplished  using  a Fortran 
language  computer  program  on  the  IBM  360/67  computer  located 
at  the  Naval  Postgraduate  School.  The  program  is  shown  in 
a separate  section  beginning  on  page  61. 
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The  calculations  involved  are  of  mean  freestream  static 
pressure,  density,  velocity  and  magnitude  of  the  oscillation 
of  mean  freestream  velocity  in  percent.  Then  the  pressure 
distribution  was  non-dimensionalized  with  respect  to  free- 
stream  dynamic  pressure.  Finally  the  blowing  coefficient 
was  calculated  for  the  upper  and  lower  blowers. 

In  addition  the  computer  program  plotted  the  pressure 
distribution  for  constant  frequency  and  blowing  coefficient 
versus  frequency  for  constant  ireestream  dynamic  pressure. 

The  pressure  distribution  plots  were  used  for  further,  manual, 
data  reduction.  The  blowing  coefficient  variation  with  fre- 
quency was  a guide  to  further1  experimentation  in  the  early 
stages,  and  for  presentation  of  results. 


V.  RESULTS 


The  results  of  this  study  are  presented  in  Figures  V.2 
through  V. 11. 

One  of  the  greatest  difficulties  encountered  during  ex- 
perimentation was  establishing  a reasonable  definition  of 
separation  in  an  oscillating  flow.  Turbulent  oscillating 
flow  does  not  lend  itself  to  a single,  unambiguous  defini- 
tion of  separation  because  the  flow  tends  to  separate  pro- 
gressively. Flow  areas  may  be  identified  that  are  attached 
always,  separated  always,  and  alternately  attached  and 
separated  during  a single  cycle  of  oscillation.  The  method 
of  visualizing  separation  with  rows  of  tufts  led  to  the 
definition  used  in  this  study.  It  was  assumed  that  separa- 
tion had  occurred  if  the  flow  reversed  itself  within  the 
boundary  layer  during  any  part  of  a cycle.  The  definition 
used  in  this  study  does  not  coincide  with  Despard's  defini- 
tion [Ref.  12]  for  separation  in  laminar  boundary  layers. 
Despard  proposed  separation  as  commencing  with  the  initial 
occurrence  of  zero  velocity  or  reverse  flow  throughout  the 
entire  cycle.  Despard's  definition  cannot  be  easily  used 
for  flow  visualization  studies.  Preliminary  testing  re- 
vealed that  below  about  1.6  cycles  per  second  the  flow  acted 
in  a quasi-steady  manner.  That  is,  the  blowing  required  to 
maintain  attachment  was  dependent  on  the  instantaneous 
velocity.  Therefore  it  was  necessary  to  provide  blowing 
equal  to  that  required  by  the  maximum  magnitude  of  the  advers 
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impossible  to  detect,  by  eye,  the  onset  of  separation  with 


the  method  in  use.  In  fact,  as  the  frequency  approached  16 
cycles  per  second  the  decision  as  to  whether  to  flow  was  or 
«as  not  separated  became  more  arbitrary  because  the  tufts 
would  oscillate  at  a frequency  close  to  t^ie  limit  of  the 
eyes  ability  to  discern  change.  For  the  reasons  cited  above 
the  study  was  continued  to  the  frequency  range  from  1.6  to 
16  cycles  per  second. 

One  other  problem  encountered  with  the  experimental  setup 
was  that  of  three  dimensional  flow  in  the  test  section.  It 
has  been  suggested  that  the  study  of  separation  on  a flat 
plate  in  a turbulent  oscillating  flow  with  strong  adverse 
pressure  gradient  could  be  accomplished  by  using  a model 
similar  to  that  used  in  this  study.  To  achieve  this  it 
would  be  necessary  to  have  the  flow  attached  across  the 
entire  plane  of  the  diverging  section.  This  should  allow? 
the  boundary  layer  growing  on  the  flat  plate  to  become  de- 
tached due  to  the  impressed  pressure  gradient.  Tufts  were 
placed,  therefore,  on  both  the  upper  and  lower  surfaces  of 
the  flat  plate,  as  well  as  along  the  diverging  plane.  Boun- 


dary layer  growth  in  the  corners  and  on  the  sidewalls  caused 
a wedge  shaped  area  of  separation  -to  appear  on  the  diverging 
plane  as  shown  in  Figure  V.l.  Tufts  close  to  tne  sidewalls 
indicated  this  corner-sidewall  separation  but  those  close  to 
the  centerline  of  the  tunnel  indicated  the  flow  was  stilll 


attached  to  the  diverging  plane.  It  was  necessary,  therefore, 
to  observe  only  tufts  in  the  center  of  the  test  section  dur- 
ing this  study.  The  flow  was  considered  attached  if  at  least 
five  contiguous  tufts  in  one  row  indicated  that  the  flow  was 
attached.  This  definition  was  somewhat  arbitrary,  but  its 
consistent  application  led  to  some  meaningful  trends . Some 
separation  did  occur  from  the  bottom  side  of  the  plate  when 
the  mean  freestream  velocity  was  138  feet  per  second  with 
high  blowing  coefficient.  Unfortunately,  this  separation 
was  not  reproducible.  The  flow  never  did  detach  from  the 
top  side  of  the  plate.  These  observations  show  that  it  will 
be  necessary  to  blow  along  the  sidewalls  to  prevent  the  cor- 
ner separation,  thereby  ensuring  more  "■  'edictable  separation 
characteristics  on  the  plate. 

To  analyze  the  results  it  was  necessary  to  determine  the 

parameters  on  which  C (the  minimum  blowing  required  to 

yR 

maintain  attached  flow)  is  dependent.  It  was  assumed  that 

the  jet  velocity  required  V.  to  attach  the  flow  was  a func- 

3R 

tion  of  the  density,  pj  mean  freestream  flow  velocity,  U^ ; 
the  length  of  the  plane  on  which  attachment  is  measured,  l j 
viscosity,  p ; the  magnitude  of  the  perturbation  of  the  free- 
stream velocity,  AU^;  the  frequency  of  the  oscillation,  w ; 
the  slot  height,  h;  and  the  effect  of  the  pressure  gradient 
as  measured  by  AP^’6,  where  6 is  the  boundary  layer  thickness 
at  the  diffuser  exit  and  AP^  the  pressure  difference  from 
the  reference  or  undisturbed  condition  to  the  exit.  This 
dependence  is  expressed  as  Equation  V-l. 
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coefficient  as 
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where  UR  is  the  ’'elocity  at  the  exit  of  the  diverging  section. 
Substituting  V-2  into  V-l  and  non-dimensionalizing  with  re- 


spect to  p,  U , and  £ leads  to 
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It  is  convenient  to  replace  y/pU^J,  by  its  reciprocal,  the 

Reynolds  number. 

C may  then  be  expressed  as 
K 
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where  b is  the  width  of  the  slot.  This  may  be  simplified  to 

V. 
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From  equation  V-3  and  V-5  the  critical  blowing  coefficient 


may  be  written  as 


. Ur  PU  £ AU  . 

CyR  = u|  3>(  ~~V~  )>(C’(  UT"  M T,  )}(V"6) 


Equation  V-6  shows  that  C is  a function  of  five  independent 

MR 

parameters.  In  the  experimental  study  the  pressure  gradient 
parameter 


s’s  jiti*#. 


<5  UE 

( r )C1  “ dj  3 

was  determined  by  the  model  and  tunnel  configuration.  Al- 
though the  pressure  gradient  was  measured  (see  later  dis- 
cussion) there  was  no  measurement  or  control  of  6 possible 
in  this  study.  This  parameter  may  be  most  important  and 

further  studies  should  be  made  while  observing  6.  Three 

5 S 

Reynolds  numbers  were  considered;  5.93x10  , 7.31x10  and 
5 

8.53x10  . These  Reynolds  numbers  corresponded  to  test  sec- 
tion freestream  velocities  of  96  fps , 118  fps,  and  138  feet 
per  second  (dynamic  pressures  of  5.0,  7.5,  and  10.0  centi- 
meters of  water),  respectively. 

Oscillation  frequency  was  varied  continuously  between 
1.6  and  16  cycles  per  second  which  resulted  in  a reduced 
frequency  range  of  approximately  .1  to  1.1.  The  percent 
perturbation  of  velocity  AU^/U^  was  a function  of  frequency, 
shutter  valve  plate  size,  and  freestream  dynamic  pressure. 

It  was  not  independently  controllable  in  this  study  since  a 
single  plate  size  was  used.  AU^/U^,  therefore  was  dependent 
on  frequency  and  dynamic  pressure.  Since  AU^/U^  did  vary  as 
the  reduced  frequency  was  changed,  the  resulting  frequency 

variation  in  C measured  could  not  be  determined  at  constant 
UR 

values  of  all  the  correlation  parameters. 

Reference  [18]  states  that  experiments  have  shown  that 

for  steady  flows  C is  not  a function  of  h/Jl  for  high  slot 

yR 

Reynolds  Numbers,  Vdi/v.  Therefore,  h/£  is  normally  omitted 
from  Equation  V-6 . Blowing  characteristics  could  vary  sig- 
nificantly with  slot  height  at  low  values  of  ti.ie  slot  Reynolds 

45 


numbers.  Slot  Reynolds  numbers  encountered  in  this  study  varied 

3 3 3 

from  6.06x10  to  11.56x10  for  the  upper  slot,  and  6.42x10  to 

11.74x10  for  the  lower  slot.  To  independently  check  the  re- 
sults of  Ref.  [18]  for  the  case  of  unsteady  flows,  this  study 

- 3 - 3 

used  two  noii-dimensional  slot  heights  of  3.15x10  and  4.2  3x10 

Figures  V-2,  3 and  4 clearly  indicate  that  frequency  has 

a definite  effect  on  blowing  requirements  to  maintain  an  at- 

5 

tached  flow.  At  a Reynolds  number  of  5.93x10  , Figure  V-2 

shows  that  there  is  no  apparent  change  in  C from  the  non- 

yR 

oscillating  condition  to  a reduced  frequency  of  0.5  where 

C increases  by  approximately  30  percent.  The  required 
yR 

flowing  gradually  decreases  through  the  remainder  of  the 

frequencies  tested  in  this  study.  For  the  Reynolds  number 

of  7.31x10^,  Figure  V-3,  the  initial  values  of  C at  a 

y R 

reduced  frequency  of  0.1  are  approximately  40  percent 
above  the  non-oscillating  condition.  At  a reduced  frequency 


of  0.13,  C drops  to  approximately  30  percent  greater.  At 

Ur 

0.48  the  blowing  required  increases  to  as  much  as  50  percent 
above  the  non-oscillating  condition  after  which  it  drops 
back  down  to  approximately  30  percent  for  the  higher  fre- 
quency oscillations.  Figure  V-4  shows  the  behavior  for 

5 ^ 

a Reynolds  number  of  8.53x10  . It  shows  that  a reduced 

frequency  of  C.08,  C is  approximately  40  percent  greater 

UR 

than  the  non-oscillating  case.  C rises  rapidly  to  100 

UR 

percent  higher  at  a reduced  frequency  of  0.12.  C stays  at 

yR 

this  level  until  a reduced  fi^cquency  of  0.29  is  reached 
where  it  drops  to  50  percent  above  the  non-oscillating 


condition.  At  a reduced  frequency  of  0.4,  the  blowing  re- 
quired then  rises  rapidly  to  70  percent  above  the  non-oscil- 
lating condition  followed  by  a gradual  decrease  through  the 
remainder  of  the  frequencies  studied. 

There  was  a significant  difference  between  C for  the 

yR 

upper  and  lower  blowers.  This  was  not  expected,  in  spite 
of  the  differences  in  slot  heights,  as  indicated  in  pre- 
vious discussion,  because  is  a momentum  coefficient. 
Therefore,  although  the  upper  blower  had  less  mass  flow, 
its  greater  jet  velocity  should  produce  the  same  attachment 
effectiveness  as  the  lower  blower.  The  upper  blower  C 

yR 

was  approximately  15  percent  higher  than  that  of  the  lower 


blower  at  R^  of  5.93x10' 


For  of  8.53x10  the  difference 


was  20  percent.  But  at  R^  equal  to  7.31x10  , the  upper 
blower  required  approximately  4 percent  less  blowing  than 
the  lower  blower.  This  inconsistency,  with  test  section 
velocity,  of  the  difference  between  upper  and  lower  C 

yR 

migh+  :r.dicatp  that  the  upper  and  lower  diverging  sections 
expe  ‘once  different  flow  environments.  Since  separation 
is  intimately  linked  with  both  previous  boundary  layer 
growth  and  disturbing  influences,  the  C differences  may 


R 


be  explained  by  the  presence  of  some  downstream  flow  asym- 
metry that  is  a func+ion  of  test  section  velocity. 

It  was  rocogn’.  ,eu  that  the  variation  in  C was  depen- 

yR 

dent  on  the  pressure  gradient  and  the  magnitude  of  the 
perturbation  of  free  '■''earn  velocity,  both  of  which  exhibit 
variations  with  fre ^ .ancy  in  this  test  setup.  The  median 


f ' 

/ 


pressure  coefficient,  , on  the  flat  plate  was  plotted  versus 


position,  x,  along  the  plate,  x/fc  equal  to  zero  corresponds 
to  the  blowing  slot  location.  Typical  pressure  distributions 
are  shown  in  Figure  V-5  for  the  cases  of  steady  flow  with  no 
blowing,  steady  flow  with  blowing,  and  oscillating  flow  with 
blowing  respectively.  The  pressure  gradient  was  determined 
by  estimating  the  slope  of  the  best  straight  line  fit  of  the 
plotted  pressures.  The  resulting  slope  was  plotted  against 
the  reduced  frequency  for  each  of  the  freestream  velocities 


run.  These  data  showed  similar  variation  of  C versus  re- 

px 


duced  frequency  and  so  are  plotted  together  in  Figure  V-6. 

An  examination  of  Figure  V-6  shows  that  the  pressure  gradient 
variation  with  frequency  was  not  the  primary  cause  of  the 


exhibited  large  variations  of  C with  frequency.  Similar 


plots  of  percent  perturbation  in  freestream  velocity,  as 


measured  by  the  hot  wire , versus  reduced  frequency  are  shown 

c 

in  Figures  V-7,  8 and  9.  The  cases  of  equal  to  5.93x10 


and  7.31x10  show  a gradual,  almost  linear,  decrease  of 


perturbation  velocity  with  frequency  but  the  of  8. 53x10s 
case  showed  the  characteristics  seen  in  Figure  V-9.  There 
was,  however,  no  indication  that  these  variations  in  the 
velocity  perturbation  were  responsible  for  t-he  shape  of  the 


large  variation  of  C with  frequency  seen  in  Figure  V-2,  3 


U 


R 


and  4 . 

The  blowing  requirement  to  maintain  attached  flow  has 
been  shown  to  be  frequency  dependent.  This  dependence  ex- 
hibits characteristics  which  suggest  resonant  behavior. 


4 8 
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This  resonance  might  be  the  result  of  flow  confinement. 

Wave  reflections*  due  to  the  tunnel  walls,  could  affect  the 
conditions  in  the  test  section.  These  reilected  waves  could 
either  reinforce  or  dissipate  the  shutter  valve  induced 
pressure  pulses.  Another  source  of  resonance  lies  in  the 
physical  makeup  of  the  tunnel.  It  was  observed  during  the 
experimental  runs  that  certain  frequency  bands  at  each  tun- 
nel test  section  velocity  made  the  wind  tunnel  vibrate  a 
great  deal,  sometimes  causing  the  test  section  to  be  travel- 
ing as  much  as  two  inches  longitudinally.  These  resonant 
frequency  bands  are  shown  in  Figures  V-2,  3 and  4.  This 
resonant  behavior  of  the  tunnel  occurred  around  a frequency 
of  approximately  8.7  cps.  Tunnel  vibrational  modes  are  seen 
to  be  influencing  the  blowing  characteristics  observed.  Dif- 
ferent test  configurations  may  exhibit  different  characteristics. 
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Figure  V.4.  Blowing  Coefficient  Required  vs.  Frequency 
Rn  = 8 . 5 3x10 $ . 
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Figure  V.7.  Velocity  Perturbation  vs.  Frequency 
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VI.  CONCLUSIONS 


From  the  results  obtained,  the  following  conclusions  may 
be  drawn: 

1.  There  is  a definite  effect  on  the  blowing  required 
to  maintain  an  attached  boundary  layer  in^a  strong  adverse 
pressure  gradient  as  oscillations  are  superimposed  on  the 
frees-!.  . . velocity  and  the  frequency  is  varied. 

2 . The  frequency  dependence  of  the  blowing  requirements 
exhibits  characteristics  which  suggest  resonant  behavior. 

3.  Considerable  three-dimensional  flow  is  produced  in 
the  present  test  setup  and  sidewall  blowing  will  be  required 
to  produce  two-dimensional  separation  on. the  plate. 
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VII.  RECOMMENDATIONS 


1.  C is  a function  of  the  five  parameters  discussed 

yR 

in  the  Results  section.  More  tests  should  be  made  in  which 
only  one  of  these  is  varied  at  a time.  As  discussed  earlier, 
measurement  of  the  boundary  layer  thickness  is  most  desirable. 

2.  The  velocity  profiles  throughout  the  flow  field  need 
to  be  studied  as  a function  of  time  while  varying  oscillation 
frequency.  To  accomplish  this  a method  for  flow  measurement 
must  be  developed  which  will  fulfill  the  following  require- 
ments*, non-interference  with  the  flow,  capable  of  measuring 
turbulent  variations,  capable  of  scanning  a relatively  large 
section  of  the  flow,  and  capable  of  resolving  rapid  changes 
in  flow  velocities  at  high  oscillation  frequencies. 

3.  An  investigation  of  the  natural  frequencies  of  the 
tunnel-model-air  flow  configuration  should  be  made  to  deter- 
mine the  effects  on  pressure  gradient,  magnitude  of  oscil- 
lating velocity  perturbation,  and  blowing  requirements  to 
maintain  an  attached  boundary  layer. 

4.  Sidewall  blowing  should  be  used  to  ensure  two  dimen- 
sional separation  of  the  flat  plate  for  further  studies. 
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COMPUTER  PROGRAM 


$ $ * * * 4>  * # $ # # $ * $ * * * % $ # $ $ 4c # ft  4 * ft  $ ft  * ft  * 4c  ft  * * ft  ft  $ $ * * ft  ft  ft  ft  * ft  ft  <i ft  ft  * ft  ft  ft 


THIS  PROGRAM  CALCULATES 


FOLLOWING 


MEAN  FREESTREAM  VELOCITY  IN  FEET  PER  SECOND 
OSCILLATION  IN  RADIANS  PER  SECOND 
FREESTREAM  VELOCITY  PERTLRBITION  IN  PERCENT 
DIMENSIONLESS  PRESSURE  DISTRIBUTION 
A.  ALSO  PLGTTED  ON  CALCCMP  PLOTTER  WITH  THE 
ABSCISSA  THE  X POSITION  IN  INCHES 
BEGINNING  AT  THE  UPSTREAM  POSITION  AND 
THE  ORDINATE  IS  THE  PRESSURE  COEFFICIENT 
JET  VELOCITY  GF  EJECTED  AIR  FOR  BOTH  UPPER 
AND  LOWER  BLOWERS  IN  FEET  PER  SECOND 
MASS  FLOW  OF  EJECTED  AIR  IN 
POUNDS  MASS /SECOND 
BLOWING  COerHCUNT  FOR  UPPER  AND 
LOWER  BLOWERS  V 

A.  PLOTTED  ON  '.ALCOMP  PLCTTER  WITH  THE 
ABSCISSA  THE  FREQUENCY  IN  RADIANS  AND 
THE  ORDINATE  THE  BLOWING  COEFFICIENT 
IS  LINE  PRINTED,  SEPARATED  BY  FREQUENCY 
LABELLED 


OUTPUT 
AND  IS 


ftftftft*ftftftftftftftftftftftftftftft  ftftftftft*ftft  ftftftftftftftftftftftftftft  ftftftftftftftftftft  ftftftftft 


INPUT  REQUIRED 

PATM  = ATMOSPHERIC  PRESSURE  (INCHES  OF  MERCURY) 

T = OUTSIDE  AIR  TEMPERATURE  (DEGREES  FAR EHEIT ) 

FREQ  * FREQUENCY  (CYCLES  PER  SECOND) 

Q = MEAN  FREESTREAM  DYNAMIC  PRESSURE  (CENTIMETERS 
OF  WATER) 

QMAX  = MAXIMUM  MAGNITUDE  OF  VELOCITY  OSCILATIONS 
(VOLTS  OUTPUT  CF  HOTWIRE  ANEMOMETER) 

QMIN  = MINIMUM  MAGNITUDE  OF  VELOCITY  OSCILATIONS 
(AS  ABOVE) 

PAMB  » AMBIENT  PRESSURE  (CENTIMETERS  CF  WATER, 

REFERENCE  HEIGHT  FROM  MULTI-TUBE  MANOMETER  BOARD 
P = PRESSURE  DISTRIBUTION,  16  VALUES  PUT  IN  1X16 
VECTOR  (CENTIMETERS  OF  WATER) 

P2U  * PRESSURE  JUST  DOWNSTREAM  OF  THE  ORIFICE  PLATE 
. IN  THE  SUPPLY  LINE  TO  THE  UPPER  BLOWER 
(CENTIMETERS  CF  WATER) 

P2L  = PRESSURE  AS  ABOVE  BUT  LOWER  BLOWER 
VFLOU  = VOLUME  FLOW  OF  BLOWING  AIR  TO  UPPER  BLOWER 
(CUBIC  FEET  PER  MINUTE) 

VFLOL  = VOLUME  FLOW  AS  ABOVE  FCR  LOWER  BLOWER 
TPU  = PLENUM  AIR  TEMPERATURE,  UPPER  BLOWER  PLENUM 
(DEGREES  FAR  ENHF I T ) 

TPL  = PLENUM  AIR  TEMPERATURE,  LOWER  BLOWER 
PPU  = PLENUM  AIR  PRESSURE,  UPPER  BLOWER 

(CENTIMERS  OF  WATER  DIFFERENCE  FROM  ATMOSPHERIC) 
PPL  a PLENUM  AIR  PRESSURE,  LCWEF  BLOWER 


INPUT  GOES  ON  FOUR  DATA  CARDS,  THE  FIRST  HAS  7 NUM 
THEY  ARE,  PATM,  T,  FREQ,  3,  QMAX,  QMIN,  PAMB. 

THE  NEXT  TWO  CARD  HAVE  THE  16  NUMBERS  OF  PRESSURE 
DISTRIBUTION.  THE  LAST  CARD  HAS,  P2U , VFLOL,  TPU, 


NUMBERS 


PPU,  P2L, 
AS  F10.4 


VFLOL , 


P2U , VFLOL,  TPU 

jMouks  are  Input 


Reproduced  (rom 
best  available  copy^ 


4 


v4T-*" 


nnnn  non  oooooooooooooooooooooooooooo 


LINE  OUTPUT 

FREQUENCY  (CYCLES  PER  SECOND) 

Q (CENTIMETERS  OF  WATER) 

VELOCITY  (MEAN  FREESTREAM  IN  FEET  PER  SECOND) 
PERCENT  (PERTURBATION  OF  VELOCITY) 

OMEGA  (OSCILATICNS  IN  RADIANS  PER  SECOND) 

CMU  (BLOWING  COEFFICIENT,  UPPER  BLOWER) 

VJU  (JET  VELOCITY,  UPPER  BLOWER  IN  FEET/SFCOND 
FLOU  (MASS  FLOW  OF  JET,  UPPER  BLOWER  IN 
POUNDS  PEP  SECOND) 

CML  (BLOWING  COEFFICIENT,  LOWER) 

VJL  (JET  VELOCITY,  LOWER) 

FLOL  (JET  MASS  FLOW,  LOWER) 

STATIONS  (PRESSURE  OISTRIBUTION^TAP  POSITIONS 
WITH  1 AS  THE  DOWNSTREAM  POSITION) 

CP  (PRESSURE  COEFFICENT  AT  STATION) 

PLOTTED  OUTPUT 

PRESSURE  COEFFICIENT  VS  X POSITION  (WITH  I BEING 
UPSTREAM) 

BLOWING  COEFFICIENT  VS  OMEGA  (PLOT  OF  CML  AND  CML 
VERSUS  OMEGA  FOR  EACH  INPUT  DECK) 


30 

35 

40 

45 

50 

55 


DIMENSI 
1PCSL( 16 
10  FORMAT 
15  FORMAT 
20  FORMAT 
1' OMEGA* 
1,T11C»  * 
25  FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
REAL*8 
REALMS 
REAL*8 
REAL*8 
READ  (5 
REALMS 
READ  (5 
DC  1000 
READ  (5 
READ  (5 
DC  100 
LPOS(K) 
PCSL(K) 
CONTI  NU 
READ  (5 


ON  CP( 16,40) ,CMU(40) , OMEGA( 40 ) , F ( 40 ) , 
),P(16) , CML (40) ,LP0S(16) , CPI ( 16 ) , QA ( 40 ) 
(7F10.4) 

(3F10.4) 

(•0»,T8,*Q*,T17, 'VELOCITY* ,T30j • PERCENT'  , 


T98 , 'CML* 


100 


,T58, 'CMU' ,T72, 'VJU* ,T84t 'FLOU 
VJL' i T123 , 'FLOL* ) 

(«0*. 10F13.5) 

(*0* ,T19, 'STATION* ,T37, *CP*/) 

('  « « I 20 , F26  « 8 ) 

('  1*  ) 

( 6A8) 

COLAST  = ',110) 

COTHE  FREQUENCY  IS  ',F10.5) 

LABEL  1/ ' CML  •/ 

LA8EL2/'CMU  •/ 

LABEL3/ ' •/ 

TITLEK12) 

,45)  TITLE1 
TITLE2(12) 

,45)  TITLE2 
J = l,  28 

, 10)  PATM ,T ,FREQ,Q,QMAX, QMIN, PAMB 
,15)  P 
K= 1 , 16 
= K 

=LPOS (K ) 

E 

,15)  P2U,VFLQU,TPU,PPU,P2L,VFL0L,TPL,PPL 


T 44  i 


CONVERT  TEMPERATURES  TO  RANKINE 

TEMP=T+459.688 
T PL=TPL  +459 • 688 
TPU=TPU+459.688 
Q A ( J ) -Q 
F(J)=FREQ 


CALCULATE  OMEGA,  MEAN 
PERCENT  PERTURBATION 


FREESTREAM  VELOCITY,  AND 


OMEGA (J )=FREQ*6. 2831 85308 
P INF= ( P ( 1 4 ) + P( 15 ) +P ( 1 6 ) ) / 3 • 0 

RHO=( <2 •0*PAMB ) -PI NF ) * PATM*. 04 120696200/ PAMB /TEMP 
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c 

c 

c 


c 

c 

c 


VEL*2.020916944D0*SCRT (Q/RHG) 

PERCCN= ( ( QMAX-QM I N ) / ( CMAX+QMI N) ) *100 

NCN  OIMENSICNALIZE  PRESSURE 

B?(i?3»UMte-p(in/Q 

200  CONTINUE 

CALCULATE  BLOWING  PARAMETERS 

RH02U* ( P2U+PATM ) * • 041206962 100/ TPU 

VLJU^9k?V?«??fo^DO.TPU»(l.O-(((2.O.PAM0)- 
1/((2.0*PAMB)-PPU) ) . 2862241256D0 ) ) 

CMU(J)  * ( FLOU*VJU  )/  ( Q*4 . 332271417  00,1 
RH02L* ( P2L +P ATM )*. 0412 06962 100/ TPL 

V JL*OSQRT ( 5454« 50724D04TPL* (l«0-(( (2.0*PAMB)' 
1/(12. 0* PA MB)“PPL))**. 2862 24 125 6D0i) 
CML(J)4fL0L*VJL)/<  0*4.332271417001 

LINE  OUTPUT 


WRITE 

WRITE 

WRITE 

WRITE 


( 6»  40) 
( 6«  55 ) 
(6*20) 
(6,25) 


F(  J) 

WKiit:  10*421  Q , V E L • PE  RC  EN  * OME  G A ( J ) * CMU  ( J ) » 
IV JUiFLOU* CML ( J) ,VJL,FLOL 
WRITE  (6,30) 

WRITE0 ( 6, 3£)6LP0S( I I , C P ( I » J) 

300  CONTINUE 
1000  CONTINUE 

; PLOTTING  ROUTINES 


NFTS=28 
EXSC= 15 • 0 
YSCL=.006 
MX=0 

C ALL  OR  AW ( NPTS  * 0MEGA,CML»1»1»LABEL1»TITLE1 
1, EXSC,YSCL,0,0*  MX, MY, 9, 15*0, LAST) 

CALLE0RAW(NPTS^MEGA,CMy,3,2,L^ 

1, EXSC ,YSCL ,0,0, MX, MY, 9, 15,0, LAST) 

WRITE  (6,50)  LAST 
DO  2000  J = l»  28 
DO  1100  1*1,16 
K= 17-1 

CPI ( I ) * CP ( K, J ) 

110°  CCNTINUEw(i6fposLt  6EL3|TITLE2t 

12.0,0.0,0,0, 1,0,8 ,8,0, LAST) 

WRITE  (6,50)  LAST 
2000  CONTINUE 
STOP 
END 


Reproduced  from 
best  available  copy. 
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